L1551NE - Discovery of a Binary Companion 
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L1551NE is a very young (class or I) low-mass protostar located close to 
the well-studied L1551 IRS5. We present here evidence, from 1.3mm contin- 
uum interferometric observations at ~1" resolution, for a binary companion to 
L1551NE. The companion, whose 1.3mm flux density is ~l/3 that of the pri- 
mary component, is located 1.43" (~230 A.U. at 160pc) to the southeast. The 
millimeterwave emission from the primary component may have been just barely 
resolved, with deconvolved size ~0.82"x0.70" (~131xll2 A.U.). The companion 
emission was unresolved (<100 A.U.). The pair is embedded within a flattened 
circum-binary envelope of size ~5.4" x 2.3" (~ 860 x 370 A.U.). The masses 
of the three components (i.e. from the cicumstellar material of the primary star 
and its companion, and the envelope) are approximately 0.044, 0.014 and 0.023 
M respectively. 



Subject headings: ISM: individual: L1551NE — stars: formation — stars: binary 
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1. Introduction 

Young T Tauri stars in Taurus have been found to have a high incidence of multiplicity, 
with the fraction of close (^lOOA.U.) companions found to be 0.40 ± 0.08 (Ghez et al. 1993; 
Leinert et al. 1993; Simon et al. 1995; Patience et al. 1998), using speckle and occultation 
techniques. Similar surveys of the Hyades cluster found a smaller fraction of 0.30 ± 0.06, 
and a still lower fraction of 0.14 ± 0.03 for G dwarfs in the solar neighborhood (Patience et 
al. 1998), suggesting an evolutionary effect and/or environmental effects during formation. 
Deeply embedded class I and class protostars cannot be surveyed using similar techniques, 
since they are not visible at optical wavelengths and are generally too deeply embedded 
even at infrared wavelengths. Millimeter /sub-millimeter- wavelength interferometry at sub- 
arcsecond resolution can resolve close companions if they have circumstellar disks. For 
example, L1551 IRS5 was recently shown to be a binary based on 7mm interferometric 
observations at the VLA (Rodriguez et al. 1998). We present here evidence that L1551 NE 
also has a binary companion. 

L1551NE (B1950 4 /l 28 m 50.5 s +18°02'10" (Draper, Warren-Smith & Scarrott 1985)) is 
a young stellar object in the L1551 molecular cloud, at a distance of 160 pc (Snell 1981). 
Discovered by Emerson et al. (1984) from IRAS data, it is the second brightest embedded 
source in the Taurus complex after L1551 IRS5, with Lw ~ 6L Q . It has a molecular outflow 
(Moriarty-Schieven, Butner & Wannier 1995). The radial density distribution of L1551NE 
has been modeled by Barsony & Chandler (1993) from 800/wn images, and by Butner et 
al. (1995) from 100/im and 200/im observations. Both found that the density distribution 
implied by the radial intensity profile was much shallower than the n(r) ~ r ~ L5 predicted 
by the "inside-out" collapse model of Terebey, Shu & Cassen (1984). Moriarty-Schieven, 
Butner & Wannier (1995) have suggested that L1551 NE may be a class source. 



2. Observations 

The observations were taken with the Owens Valley Radio Observatory millimeter-wave 
array 1994 November 12 and 1997 January 17, at a wavelength of 1.3mm (230.65 GHz). The 
1994 observations were taken with the array (consisting at the time of six 10m antennae) 
in the compact "Z" configuration (maximum baseline 115m, minimum baseline 30m), and in 
1997 the data were taken using the extended u h" configuration (maximum baseline 241.7m, 
minimum baseline 35m). The effective resolution of the synthesized beam was ~2" and ~1" 
for the two configurations. 
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The observations were made in "snapshot" mode, with a series of ten minute integrations 
interspersed with observations of other sources and the phase calibrators (0528+134 and 
3C84), to yield a total integration time of ~1 hour in the "/" configuration and ~2 hours 
in the u h" configuration. The system temperatures were ~1200K and ~1600K respectively. 
Neptune was used as the primary flux calibrator. The continuum bandwidth is 1GHz. 

The data were reduced using the standard reduction package "mma" , and then exported 
to AIPS which was used to generate and clean images (using the 'IMAGR' task and uniform 
weighting). The clean box included only a small area covering the main source and the 
apparent binary companion. The angular resolution of the cleaned image was 1.29" x 1.07" 
at position angle -74°. 



3. Data 

Figure 1 presents a contour image of L1551NE. Clearly seen are the primary source 
(henceforth source A) at (B1950) 04 fe 28 m 50.552 s 18°02'09.85", and another weaker source 
1.43" to the southeast (source B). Emission also appears to extend to the north and to the 
east of the primary source, and surrounding source B. 

To verify that we are seeing real sources and not phase instabilities, we generated images 
using the same techniques of another source, IRAS 04169+2702, which was observed during 
the same two days and interspersed with L1551NE. Phase errors would be manifested as 
"anomalous" sources or structures. No such anomalous sources are seen in the image of 
IRAS 04169+2702. In addition, we cleaned the "/" configuration data separately from the 
"/i" configuration data. The high-resolution data clearly had two peaks, while the low- 
resolution data showed an extended disk-like structure with long axis through the line joining 
the two sources. Thus we believe that the structure seen here in L1551NE is real. 

In Table 1 we present the positions, sizes and flux densities of the sources. The single- 
dish flux density at 1.35mm is also shown. 



4. Discussion 

There are three distinct components apparent in the image shown in Figure 1; a brighter, 
possibly extended source at the field center (source A), a weaker, probably point-like source 
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(B) 1.43" south-east of A, and diffuse, low-level emission which appears to surround both 
sources and extend ~2" to the northwest and east of A. 

Sources A and B were fit with elliptical gaussians using the AIPS task JMFIT. The 
primary source A was found to have a size 1.53"xl.28" (i.e. it may have been slightly 
resolved with a size ~ 2a larger than the beamsize), peak intensity ~ 0.33 Jy and integrated 
intensity ~ 0.47 Jy (i.e. ~40% (~ 3a) larger than the peak intensity). If it is slightly resolved 
(deconvolved size ~ 0.82"x0.70"), then its size is ~131xll2 A.U. at a distance of 160pc. 
Source B is located 1.43" (229 A.U. at 160pc) south-east of A, and its size and intensity are 
consistent with it being unresolved, i.e. <100 A.U. We estimate the mass within each source 
using M D = J^fj f^-, and display these in Table 1. We used the integrated flux densities, 
assumed a distance of 160 pc (Snell 1981), used dust temperature of 42K (Moriarty-Schieven 
et al. 1994), and dust emissivities from Hildebrand (1983) = 0.1(A/250) /3 g cm" 2 , 
assuming (3 = 1 (Moriarty-Schieven et al. 1994), and have assumed a gas to dust mass ratio 
of 100. The derived masses are approximately 0.044 and 0.014 M for the emission from 
source A and B respectively, and O.O22M for the envelope. The total mass for the two 
sources plus envelope is ~0.078 M Q . 

Rodriguez, Anglada & Raga (1995) imaged the L1551NE region at A3. 5c m and found a 
continuum source located within 1.5er of our source A. Rodriguez et al attribute the A3. 5cm 
emission to shocks associated with the outflow (Moriarty-Schieven, Butner & Wannier 1995). 
Thus source A is the likely origin of the L1551NE outflow. Rodrguez et al. found a possible 
second source located ~0.6" east of A, which is not coincident with our source B, for which no 
A3. 5cm emission was detected. If their source 0.6" east of A is real, this would suggest that 
L1551NE is at least a triple star system. However, another protostar capable of generating 
an outflow jet should have had a circumstellar disk large enough to be detected by our 
observations. Possibly the eastern A3. 5cm source is due to a jet from A but offset from the 
source, or is a background object. 

The single-dish flux density (Butner et al. 2000) found using the 14m JCMT telescope 
(FWHM 20")) is not significantly different from the total intensity found in our image of 
L1551NE (primary beam FWHM 28"). Thus only a small amount, if any, of the single- 
dish flux can have been resolved out by the interferometer. However, considerable extended 
emission was seen at ~850 fim by Barsony & Chandler (1993) and Moriarty-Schieven et al. 
(1999), and at 200 /im by Butner et al. (1995). Barsony & Chandler and Butner et al. 
modeled this extended emission as an envelope whose radial density distribution decreases 
very slowly with distance from the protostar. 

Weak, low-level emission can be seen extending ~l-2" to the north and east of source 
A, and perhaps encompassing source B. This extended structure has a disk-like appearance, 
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of dimension ~ 5" x 2" (~800x300pc with long axis at position angle ~-12°. This is roughly 
perpendicular to the axis of the conical reflection nebula emanating from L1551NE (Draper 
et al. 1985; Hodapp 1995), and hence of the molecular outflow (Moriarty-Schieven et al. 
1995). This disk-like structure may represent a circum-binary disk. 

JAP was supported by a Hawai'i Space Grant College Fellowship which is funded by the 
NASA Undergraduate Space Grant Fellowship program. The Owens Valley millimeter-wave 
array is supported by NSF grant AST-96-13717. 
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Table 1: Source Parameters 



Source A 




Peak Intensity 


0.333 ± 0.022 Jy/beam 


Integrated Intensity 


0.473 ± 0.049 Jy 


Position (B1950) 


04 A 28 m 50.559 s 18°02'09.91" (± 0.1") 


Size 


1.53" x 1.28" ± 0.1" P.A. 104° ± 15° 


"Deconvolved" Size 


0.82" x 0.70" (131 x 112 A.U. at 160pc) 


Mass" 


0.044 M 


Source B 




Peak Intensity 


0.146 ± 0.022 Jy/beam 


Integrated Intensity 


0.197 ± 0.048 Jy 


Position (B1950) 


04 /l 28 m 50.604 s 18°02'08.64" (± 0.17") 


Mass" 


0.014 M 


Circumbinary Disk 


Integrated Intensity 


0.233 ± 0.048 Jy 


Size 


5.51" (± 0.44") x 2.59" (± 0.22") P.A. 2° (± 14°) 


Deconvolved Size 


5.4" (± 0.3") x 2.3" (± 0.2") P.A. 2° (± 10°) 


Mass" 


0.022 M 


Total Integrated Intensity 


Integ. Intensity 


0.851 ± 0.084 Jy 


Mass" 


0.079 M 



Single Dish Intensity 6 

19" FWHM beam 0.83 ± 0.03 Jy 
Mass" 0.078 M fi 



"Assuming T d =42K, M 9 /M d =100, D=160pc. 

6 From Butner et al. (2000). Obtained with 14m JCMT (20" FWHM)). 
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Fig. 1. — Contour image of L1551NE. Source A is the brighter source at the field center, B 
is 1.43" to the southeast. The two lowest positive contours are at 0.049 and 0.07 Jy/beam, 
and increment at intervals of 0.035 Jy/beam. The negative contours are -0.049 an -0.07 
Jy/beam. RMS noise of the image is 0.023 Jy/beam. The synthesized beam (FWHM) is 
shown in a box in the lower-right corner of the image. 



